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Abstract 



The charge-dependent realistic nuclear Hamiltonian for a nucleus, composed 
of neutrons and protons, can be successfully approximated by a charge- 
independent one. The parameters of such a Hamiltonian, i.e., the nucleon 
mass and the NN potential, depend upon the mass number A, charge Z and 
isospin quantum number T of state of the studied nucleus. 
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I. INTRODUCTION 



The model of the atomic nucleus as a system composed of neutrons and protons has 
been studied for over sixty years and has been found to provide a more-or-less adequate 
description. By considering the constituents of the nucleus as structureless identical fermions 
(i.e., nucleons), we obtain a very useful formalism, which leads to significant simplifications 
of our understanding and description of nuclear structure. These simplifications follow from 
the fact that the neutron and the proton have similar masses and that the potentials of the 
strong neutron-neutron (nn), neutron-proton (np) and proton-proton (pp) interactions are 
also very similar. Over the intervening years, simple models without correlations, such as 
the nuclear shell model, have been unable to distinguish between different parametrizations 
of the nucleon-nucleon (NN) interaction or to detect details in the definition of the mass of 
the nucleon. 

During the last decade, however, theoretical techniques for calculating few-nucleon sys- 
tems and the lightest atomic nuclei have advanced to the point that the quality of the 
realistic NN potential is starting to play an important role |IJ . This necessarily involves the 
introduction of charge dependence into the NN interaction. The most-recent realistic poten- 
tials, such as the Nijmegen [Q, Argonne [[| or CD-Bonn are charge-dependent. On the 
one hand, the charge dependence of the nuclear Hamiltonian complicates the description, 
but on the other, the deviations of charge- dependent potentials from the charge-independent 
ones, although non-negligible, are not very large. As a consequence, modern calculations of 
light nuclei use as a starting point a Hamiltonian preserving the isospin quantum number 
and include charge-dependent effects perturbatively later on. The charge dependence of 
the kinetic and potential energies of a real nucleus, however, can only be successfully taken 
into account perturbatively, when the zeroth-order, charge- independent, isoscalar nuclear 
Hamiltonian is as close as possible to the charge-dependent one. We present here a method 
of construction for such a Hamiltonian. 

In Sections II and III we show that the kinetic-energy operator of an atomic nucleus 
composed of neutrons and protons can be precisely approximated by a corresponding opera- 
tor for a nucleus composed of nucleons, when the nucleon mass is defined as the mean value 
of the quantum-mechanical operator for the nucleon mass. The nucleon mass obtained in 
this way depends on the mass number A and charge Z of the studied nucleus. 

In Section IV we investigate the problem concerning the definition of the NN interaction. 
For the sake of simplicity, we consider potentials, defined separately in different two-nucleon 
channels, such as, the Reid'68 |J and Nijmegen || potentials. The generalization of results 
given by other realistic charge-dependent potentials, such as the CD-Bonn Q and Argonne- 
18 |J potentials, is straightforward. It is shown that the best NN potential is a linear 
combination of the nn, np, and pp potentials with coefficients, dependent not only on the 
nucleus studied, i.e., A and Z, but also on the isospin quantum number T of the state 
under consideration. So, the isoscalar Hamiltonians that we introduce are different for 
different nuclei and even for different states of the same nucleus. The derived values for the 
nucleon mass and the NN interaction must satisfy minimal requirements (e.g., for systems 
composed only of neutrons, the nucleon mass formula must give the neutron mass and the 
NN interaction must reduce to the potential of the neutron-neutron interaction) and must 
also be equal to the corresponding values for non-trivial cases suggested for three-nucleon 
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nuclei in earlier publications @, @. 

In Section V we present numerical tests of our recommendations and give our conclusions 
in Section VI. 



II. THE MASS OF THE NUCLEON 



The nonrelativistic kinetic energy operator of the atomic nucleus equals 



fc2 ^ p.2 A »-2 

F i=l i=Z+l 

where m p and m n are the proton and the neutron masses, respectively; M = Zm p + Nm r , 
is the mass of the nucleus and 



R = li fe m ^+ E 



.4 

m n Yi 



,i=l i=Z+l 

is the radius-vector of the center-of-mass of the nucleus . 

The kinetic energy operator of the nucleus, as a system of nucleons is simpler: 

*° = -|;X>? + ^, (2) 

where m is the nucleon mass and 

1 A 

i=l 

is the center-of-mass radius-vector of the nucleus composed of nucleons. 

Obviously, the operators K and K° are different. In particular, the expression for K° 
contains the parameter m, which, in many applications, is defined as an average of the 
proton and neutron masses; i.e., 

m n + m p 

m = . (3) 

On the other hand, calculations for a neutron and proton system involve a transformation 
to relative and center-of-mass coordinates and a mass equal to twice the reduced mass 

2m n m p 
m n + m p 

While in the three nucleon case, it is sometimes assumed that 3 H and 3 He are composed of 
identical particles, each of mass | (2m n + m p ) , | (m n + 2m p ) [fj , respectively. 

At first glance, it looks as though we can define the optimal kinetic energy operator, 
which is as close as possible to a charge-dependent one, yet charge- independent, by properly 



3 



defining the mass of the nucleon in K°. However, the nucleon mass is not really a free 
parameter. It has to be defined as the mean value of a nucleon mass operator: 



l a l A 

m = -^m(i) = -^ {m n P n (i) + m p P p (i)) , 



i=i i=i 
where to (i) is the mass of an i-th nucleon and 



(5) 



are projectors upon proton and neutron states. 

The mass of the nucleon in the nucleus (A, Z) can be defined exactly, because in the 
A-nucleon state, with a defined isospin projection, M T = | (N — Z) : 



m = (M T \m\ M T ) 



M' 



T 



1 A 

i=l 



M 7 



m n + m p (m n - m F 



2 

+ 
+ 



E 



M 7 



i=l 



2 A 

m n + m p (m n - m p ) 



M T , 



or, in other words, 



m = — (Nm n + Zm p ) . 



(6) 



This definition of the nucleon mass coincides with that used in the above-mentioned three- 
nucleon case. For light nuclei with N = Z , it equals the mean value of the neutron and 
proton masses. For two nucleon systems it coincides with twice the reduced mass in cases 
of two protons or two neutrons and is very close to twice the neutron-proton reduced mass 
because 



pi — m 



2m n m p 



m n + m p 2 4 MeV 
1 = - m S z ?s -4.5 x 10 — 



where the widely used small parameter is: 

?7i „ — m. 



m n + m p 



^ w 6.887 x 10" 4 . 



(7) 



Obviously, the difference between \i and m is negligible. Realistic NN potentials, defined by 
using a reduced mass for the two nucleons, cannot change significantly when the nucleon 
mass to, defined for the two-nucleon system, is used instead of ji. 
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III. INTRINSIC VARIABLES AND OPERATORS 



A direct comparison of the operators K° and K is impossible due to different neutron 
and proton masses and different definitions of center-of-mass operators. On the one hand, 
both operators, K° and K, have redundant variables, because they are functions not only 
of one nucleon but also of center-of-mass variables. On the other hand, both operators 
are intrinsic; this means dependent only on intrinsic, translationally invariant variables and 
independent of center-of-mass radius vectors. In order to compare them, it is necessary to 
introduce intrinsic variables. 

The Jacobi variables r] , rji, ...r/A-i for particles with equal masses are well-known (see 
H , || and references therein). They can be defined by a corresponding Jacobi tree as an 
orthonormal set. The expressions for rj a in terms of and vice versa are as follows: 

A A-1 
i=l a=0 

Here a is an orthogonal (A x A) matrix : 

A A-1 

«a,ifl/3,i = 5 a ^, Cba^aJ = (9) 

i=l a=0 

The translational invariance of the intrinsic variables rji, ...T]a-i along with orthogonality 
gives the following conditions for the matrix a : 
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VAS a ,o, ao,i = -!=■ (10) 



Using these variables, we can rewrite the operator O) as follows: 



fc2 A-i 

a=l 

The analogous presentation for K requires Jacobi variables for particles with different 
proton and neutron masses. Let us introduce Jacobi variables in the general case, in which 
the masses of all particles are different. The reduction to protons and neutrons will then be 
straightforward. These variables can be defined by the same kind of Jacobi tree, as in the 
previous case. 

The Jacobi tree by definition has (2 A — 1) vertices, A of which are vertices of the first 
degree (the degree of a vertex is defined as the number of edges matching in this point). They 
should be arranged in a line and marked with the one-particle radius- vectors r 1; r 2 , i\4 
. In case of particles with different masses, however, the vertices of the first degree of the 
tree have to be marked not only by one-particle radius-vectors but also by corresponding 
masses mi, m 2 , ttia- The remaining (A — 1) vertices (situated below the first group) 
determine the Jacobi coordinates £ a (a = 1,2, ...,A — 1) .The degree of each equals three, 
except for the deepest one which equals two. Contrary to the case of equivalent particles, the 
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orthogonal transformation can now be defined only between Jacobi variables £o>6i> ••• ) £a-i 
and modified one-particle variables x 1; x 2 , x A , where Xj = ^/m^ri. For the a-th vertex the 
Jacobi coordinate is 



\ 



(«)„» 

R 



m T m 



ni 



(«) 



+ m 



(a) 
R 



III 



(a) 



E 

ie{L} 



m,Xi — 



(a) 
R 



E 

ie{R} 



m,;X,; 



1,2,..., A 



(12) 



Here {L} is a manifold of the first-degree vertices, which could be reached while moving 
from the a-th vertex upwards along the left edge, while {R} denotes the same for the right 
edge. Obviously, the new Jacobi variables, £ a (a = 1,2, A — 1) , like the corresponding 
old ones, r] a (a = 1, 2, A — 1) , are translationally invariant, because they are defined as 
properly normalized radius vectors between the centers-of-mass of two subtrees, defined by 
left and right edges of the tree, as they are seen from the a-th vertex. The masses of left 
and right Jacobi clusters are defined in the following way: 



in 



(a) 



E 

i£{L} 



rrii 



m 



(a) 



E 

i£{R} 



rrii 



This set has to be completed by a zero Jacobi variable, situated on the left side of a second- 
degree vertex and proportional to a center-of-mass radius vector: 



1 A 



(13) 



where M = Ylt=i m i * s the mass of the system. In general, 



= ^ K,i (m x , m 2 , . . . , m A ) x< 



(a = 0,1,...,A-1) 



i=l 



(14) 



where b is an orthogonal (A x A) matrix. In the case of equivalent particles, this transforma- 
tion coincides with that defined in Eq. (|8|) multiplied by y/m. The proof of the orthogonality 
of these transformations is given in Appendix A. 

The tree sample for five particles is presented in Fig. 1. The matrix b for this tree equals 



/ 



V 



i mi 
M 



mim4 5 
Mm-123 

mi2«il23 



m 2 
mi2 



M 



"12 "145 



Mm-i 



m2m3 
mi2«ii23 



mi 
mi2 



M 



Mm-i 



mi2 
"1123 



M 



"M?"123 
Mr"45 









"145 



r mr L 
M 



\ 



m 5 ini23 
M"l45 








"14 
m 45 



where mj,„ k = Y^i=j m i\ M = m i2345- 
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Due to the orthogonality of the transformation from one-particle variables Xj., X2, x^ 
to a system of Jacobi variables £o>6i> ...,£4-1 , the transformation from one set to another 
set of Jacobi variables, defined by different trees, could be determined in the same way. 

The kinetic energy operator ([!]) in terms of these new variables can be rewritten as 



K 



h 2 A ~ L 



(15) 



a=l 



A comparison of operators Eq. flTTD and Eq.([l5|) requires the transformation matrix between 
the sets £ 0; 61, ■••> £a-i an d r} , rji, ...t)a-i ■ The result cannot depend on the choice of the 
Jacobi tree, so we can choose the simplest one, given in Fig. 2, where protons are marked 
by the numbers 1, Z and neutrons by the numbers Z + I, Z + N = A . In such a case 



[m 1 ,m 2< 



m A ) Xj = ^ h ^ ( m i> m 2 ,...,m A ) 



»=i 

A-l 



A-l 
(3=0 



^ ( ^b a ,i (m 1 ,m2,...,m A ) 

13=0 



midpj rip. 



Due to the multipliers 



the matrix of the last transformation is not an orthogonal 



matrix. This gives the following expression for K in terms of variables r\ 



K 



2 

M 



ZN(m„—m p y 



V 2 + 

v Jin 1 



2\IZN{m n -m p ) 



AMm n m p ' m 
72 + _1 S^Z 



Am n m p 



(V 



-V 2 + — V V 2 +— V 

Am n m p Vi m p i-^a=2 v n a 1 m n ^ 



A-l 



V 2 

Z+l %a 



So, the neutron and proton mass difference causes the intrinsic kinetic energy operator 
to depend upon the center-of-mass variable. The center-of-mass starts to move in an atomic 
nucleus composed of nucleons instead of real protons and neutrons. However, this effect 
is not significant, because the first term, proportional to V 2 , has a very small coefficient, 
while the expectation value of the second term, proportional to (V^ ■ V^J , vanishes either 
when the excitations of the center-of-mass of a nucleus are under control or when the wave 
function is intrinsic. 

The difference in the kinetic-energy operators for neutron-proton and nucleon nuclei is, 
thus, given by 



2m 



+ 



K-K° = 

Mm _ ^ \ y 2 



^ZJVKjf^f y2 _|_ 2$^- ZN m 



Mmn'iriv 



+ 



-IE 



V 2 

'a=2 v Vo 



m n m p 



+ ( jul - A y^- 1 v 2 

^ y mn L J l^ a =Z+l V r ]a 



where 5 is the parameter (^). 

From the above expression, it is obvious that the mean value of the difference in kinetic 
energies, estimated using intrinsic wave functions, equals the sum of the mean values of the 
last three terms and can be rewritten as: 



{SJUM T I K - K° | SJUMt) 



Mm 
Am n m p 



x (EJTIMt 



m 
m 



v 



l){Z-l) 



m 

m. 



n 



-l)(N 



& v.2 

V 

2m v ° 



SJUMt 



(16) 
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because the mean value of is independent of a (see Appendix B). The expression given 
in square brackets equals 

ZN (A - 1) (m n - m p f _ AZN (A - 1) a 
A 2 m n m p A 2 

while the matrix element has the same order as the kinetic energy of a single nucleon in a 
light nucleus ( < 40MeV, see fL(| ). Consequently, 



(K - K°) < 10~ 5 MeV. 

This means that by using the defined value of the nucleon mass (H) , we can account for the 
n-p mass difference in the expression for the intrinsic kinetic energy without any problems. 
In cases, when the mean value of the neutron and proton masses (0) is used instead of the 
defined value of the nucleon mass, this difference of expectation values (|1(J) is proportional 
to 5 and can cause a noticeable difference in the kinetic energies K and K°. 

IV. THE POTENTIAL OF THE NUCLEON - NUCLEON INTERACTION 

The Hamiltonian of an atomic nucleus with a charge independence and charge symmetry 
breaking NN potential has the form 

H = K + V, 

where K is the translationally invariant operator of kinetic energy, Eq. ([!]), and 

A 

y= i y pp (*> fc ) p p w p p w + Vn p (*> fc ) Pn w p p ^ + v ™ (*> fc ) p « ( f ) p « 

i,fe=l(i<fe) 

is the potential-energy operator (the Coulomb interaction is included in V pp (i, k)). 

It is useful to employ the isospin formalism and the partial-waves expansion in order to 
simplify the definition of the NN potential, i.e., 

V (i, k) = V (r ik (TiTi<JkTk) 

= Yl VJntmt ( r * fc ) P ^ ( e ik<PikViTia k T k ) , (17) 

j-ntmt 

where Yi k = {rik6ik<Pik) = r « — r fc is the difference of radius-vectors of the i-th and the k-th 
nucleons; are the spin and isospin variables of the i-th nucleon, respectively; jirb are 
the quantum numbers - total momentum, parity and isospin - of the two nucleon state; and 
Pjirtmt {9ik4>ik&i T iC'k' T k) is the projection operator on this state. The correspondence between 
our notation, j7rt, and conventional spectroscopic identifiers 2s+1 Lj is: 

0+1 - 1 S , (T1~ 3 P , l + ~ 3 St - 3 £>i, n~ 3 Pi, l^O - 1 Pi,.... 

Also, m t stands for the two-nucleon isospin projection; m t = 1 corresponds to two neutrons; 
m t = 0, to neutron-proton; and m t = —1, to two protons. Traditionally, NN potentials 



S 



are constructed by fitting np {m t = 0) data for t = states and either np {m t = 0) or pp 
{m t = — 1) data for t = 1 states fKj . 



In general, the realistic potentials can be different for m t = 0, ±1 in each two- 
nucleon channel with t = 1. To some extent it holds even for charge-symmetric poten- 
tials of the strong interaction, because the Coulomb interaction must to be included in 
yjn(t=i)(mt=-±) ( rifc ), i n channels with ix = ( — l) J the potential V^ tmt {r%k) is one function 
of , when t = 0, and three different functions of r^, when t = 1 (the best example being 
1 So). In channels with tt = (— 1) 3+1 and t — 0, the potential V^ ntmt (r ik ) can be given as three 
different functions of r ik , namely entries in a second-order symmetric matrix (such as the 
parametrization in the channel l + ~ 3 Si — 3 Di , as given for the Reid68 H and Nijmegen 
potentials) or nine different functions of r^, corresponding to different values of m t when 
it = (— 1) J ' + and t — 1. The first channel of this kind is 2~1. This parametrization, however, 
is performed as three different functions of r^. As mentioned earlier, they obviously are 
different for m t = — 1 in comparison to m t — or m t = +1, due to Coulomb interaction of 
protons. 

The charge-independent Hamiltonian of an atomic nucleus composed of nucleons is 

H° = K° + V° (18) 

with the kinetic energy operator, defined in Eq.(|2|), and the charge-independent potential- 
energy operator 

v°= fl v °^ k )> ( 19 ) 

i,k=l(i<k) 

whose partial-wave expansion is given by 

V° (i, k) = V° (rikCTiTiCTkTk) 

= ^ V jnt (r ik ) P j7Tt (Oik^ikCrinakTk) . (20) 

jnt 

As was shown above, there are no problems with regard to the charge- dependence of 
the kinetic energy, because by taking the proper value of the nucleon mass, one can use 
instead of K in the expression for H, due to the negligible difference between these two 
kinetic energy operators. In such a case, the Hamiltonian H can be expressed in the form 

A 

h= y, h ^ k ^ ( 21 ) 

i,k=l(i<k) 



where 



ft 2 2 

h(i, k) = -— — (Vi - Vfc) +V(r i -r k ,a i T i a k T k ). (22) 



The exact expression for an arbitrary eigenvalue of such an operator is ( see [11], [12 
and references therein) : 
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£ jum t = {H) JUM T = £H£i ij £ ^Q £)i7rtmt (5 JHM T ) , (23) 

£,j7rimt 

where the sum runs over all the two nucleon states jirt and all eigenvalues e of the 
Reduced Hamiltonian (RH) operator 

h = h (A - 1, A) = — ^A r + V (r, (Ta-iTa-i^t^) (24) 
with r = r A ~\ — r^. 

The Qej-Ktmt {SJUMt) are diagonal entries of the intrinsic density matrix in terms of the 
RH eigenfunctions: 

Qejntmt {£ J~RM T ) = Qs,jirtmt,£,jirtmt {£ J^M T ) , (25) 



Qe,jntm t ,e',j'n>t'm' t (£JRM T ) = ^ J d£l...cZ£^_ 2 c££^_ 1 G?£^_ 1 

aiTi...a A _ 1 T A _ 1 a A T A a' A _ 1 T , A _ 1 a' A T , A 

X ^t m (U-lVA-lTA-lVATA) 

X^sjilMt (€i~-€a-2£a-i<tiTi...o A-lTA-lO at a) 
X ^£JUM T {^■■■^A-2C A ~i (r i T i--cr' A ^ 1 T' A _ 1 cr' A T' A ) 

*t Wml (W-iVa^) > (26) 

where jnA/ T (Ci---Cyi-i (7 i' r i--- cr ^TA) is the intrinsic wave-function of the nucleus. 
According to the definition of the density matrix 

Qe,jntm t (£JUM T ) > 0, and Qejmmt (£ JHM T ) = 1. (27) 

e ,jirtmt 

These are probabilities of definite states of two nucleon relative motion, ejntm t , in the 
state under investigation, £JUMt, i.e., the state specified by the exact quantum numbers 
- energy, total momentum, parity and isospin projection of the nucleus (A, Z) . 

The expression for the eigenvalue of the Hamiltonian 

A 

H°= h °^ k ^ ( 28 ) 

i,k=l(i<k) 

where 

h 2 

h° (i, k) = -—— (V, - V k f + V° (Ti - r k , aiT i( j k T k ) (29) 
2mA 

is very similar to that given by Eq. (f23|): 

g mT = {H oynT = A(A-l) ^ £jntQ ^ {SJUT) (3Q) 

E,jirt 
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but is slightly simpler, because the corresponding NN potential is charge-independent. As 
a consequence, the total isospin T is an exact quantum number. Moreover, in such a case, 
the eigenvalues are the same for all nuclei of the isospin multiplet and, hence, independent 
on Mt ■ Here 

h o = h o {A _ 1; A) = -—A r + V° (r, a A ^T A _ x a A r A ) . (31) 



A straightforward comparison of the eigenvalues, Eqs. ( p3|) and (0), is impossible, 
as eigenfunctions of the Hamiltonians, H and H°, do not coincide. Eigenvalues and 
eigenfunctions of the corresponding Reduced Hamiltonians (Eqs. (|4j) and (|3l|)) are also 
different. To compare them, it is necessary to define kind of equivalent potential V°, for 
which we are looking. From partial-wave decompositions of the potentials, Eqs. (|17]) and 
(POD, it follows, that the realistic potential in the states with t — is the same in both 
cases, because m t can take only the value zero. If t = 1, there are three different pairs of 
nucleons (m t = 0, ±1) and the NN potential V^ tmt (r^) is different in each case. Therefore, 
the symmetric, charge-independent potential in an arbitrary two-nucleon state can be defined 
as the sum of these charge-dependent potentials: 

V J7Tt (r ik ) = <hrn t V j * tmt (r ik ) , (32) 

where the c tmt are normalized coefficients: 

^2c tmt = l. (33) 

mt 

The optimal values of these coefficients are (see Appendix C): 

coo = 1, (34) 
{A - 2M T ) {A - 2M T - 2) = 4Z (Z - 1) 

Cl_1 3A(A-2)+AT(T + l) 3A (A - 2) + AT (T + 1) ' 1 j 

A (A - 2) + AT (T + 1) - 8M T _ 8NZ - A (A + 2) + 4T (T + 1) 
Cl ° ~ 3A (A -2)+ AT (T+ 1) ~ 3A (A - 2) + AT (T + 1) ' ^ 
(A + 2M T ) {A + 2M T - 2) = AN (N — 1) 

~ 3A (A - 2) + 4T (T + 1) ~ 3A (A - 2) + AT (T + 1) ' ^ 

These values, obtained as a result of consideration of probability distributions, have a 
rather simple interpretation. Having in mind the definition of marginal probabilities, as 
given in fTT||, and the definition of the c tmt , as given in (C5), we find that these coefficients 



equal the number of pairs of nucleons with the isospin and corresponding projection quantum 
numbers, t and m t , divided by the number of pairs with the isospin quantum number t. 
Finally, the charge independent Hamiltonian can be expressed in the form 

A 

H°(A,T,M T )= h°(A,T,M T ;t,k), (38) 

i,k=l(i<k) 

where 
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h° (A, T, M T ; z, fc) = - , (V, - V,) 2 + F° (A, T, M T ; r, - r fc , a^n) . (39) 

2m (A, Mt) A 

The eigenvalue of H° (A, T, M T ) is very similar to that given by Eq. (|30|) , but now has 
the form 

S mT (A, T, M T ) = (H° (A, T, M T )) J £ J = - {A ~ l) ^ (A, T, M T ) Q £Jnt (SJUT) . 

e,jirt 

(40) 

From the last expression we observe that our result for a charge-independent Hamilto- 
nian, which includes charge dependent effects (e.g., Coulomb interaction, charge dependence 
of the potential and the neutron and proton mass difference), is a function of the parame- 
ters A, T, and Mt- Comparing with the exact expression (|23|) , we observe that isospin now 
appears as an exact quantum number and the density matrix is independent of Mt- Conse- 
quently, the density matrix is the same for all members of an isospin multiplet, as it must 
be for a charge independent Hamiltonian. Contrary to the commonly used expression for 
charge- independent potentials, i.e., Eq. (|3~0"D, the present effective Hamiltonians (Eq. (|38|)) 
are different for different nuclei of an isospin multiplet and for states of the same nucleus 
with different values of total isospin T. 



V. NUMERICAL TESTS 

The effective Hamiltonians (Eq. fl3"8p), obtained above, agree with the well-known, trivial 
isoscalar Hamiltonians for nuclear matter and for states of finite nuclei with total isospin 
T = , because, in this case, the nucleon mass equals half of neutron and proton masses and 
the effective potential is the average of the three potentials. However, even in this case, the 
recommendations given above allow us to express in isoscalar form the Coulomb interaction 
of the protons. In the light-nuclei region for isospin multiplets with T ^ 0, the values 
of the nucleon mass and the effective potential are different for nuclei with different Mt 
values. 

For the triton (A = 3, T = 1/2, Mt = 1/2), for example, the mass of the nucleon is 

2 1 

m = -m n + -mp , 

and the equivalent charge-independent potential is 

i o 

yin*=l) ( r ) — tyMt=l)(rrH=0) ( r ) _ yjir(t=i)(mt=l) ^ 

3 3 

or simply 

V{r)= l -V'; 1 {r) + 2 -V nn {r) , 

which agrees with the results obtained in |J using the same nucleon mass and in for the 
charge- independent potential in the 1 S state. 
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For 3 He (A = 3, T = 1/2, Mt = —1/2) the mass of nucleoli equals 



1 2 
m = -m n + -m p . 

The equivalent charge-independent potential is 

i p 

yjir(t=l) ^ _ ±yjir(t=l)(m t =0) ^\ ± yj7r(t=l)(m t =-l) ^ 

3 3 

with the Coulomb interaction included in the last term. 

In order to test the results obtained in section IV, we performed several calculations 
in the framework of the large-basis, no-core shell- model approach [TjJ, using interactions 



that directly break isospin as well as the corresponding NN interactions given by formulas 
([34l)- (|3~7D . The effective interactions were derived from the free np, pp, and nn interactions. 
These effective interactions then served as input for the large-basis shell-model calculations 
of the A = 3 system, e.g., 3 H and 3 He, and of the A = 6 system, e.g., 6 Li, 6 He, 6 Be. 
In the present calculations we employed the CD-Bonn |4j] potential, which includes isospin 
symmetry breaking terms. In addition, the Coulomb potential was added to the proton- 
proton interaction. To get energy convergence, one needs iV max > 30 for the three-nucleon 
problem. However, for the purpose of the present work it is enough to perform the shell- 
model calculations in a sufficiently large, but still easily accessible model space. We chose 
iV max = 8, which corresponds to 8hQ excitations above the unperturbed ground state for 
the A = 3 system and QMl excitations for the A = 6 system. For the harmonic-oscillator 
frequency, we picked Ml = 19 MeV for the A = 3 system and Ml =17.2 MeV for the A = 6 

1 2 

system. These choices follow from the phenomenological formula Ml = 45A~s — 25A~3 
MeV. 

The wave-functions obtained in the no-core shell-model approach satisfy all the require- 
ments necessary for realistic nuclear wave functions, i.e., the excitation of the center-of-mass 
of the nucleus is under control and they correspond to exact eigenvalues with good quan- 
tum numbers J, II, (and T in the case of the isospin invariant effective interaction). The 
results of our shell-model calculations are presented in Table I. These results demonstrate 
that, in general, the suggested definitions for the nucleon mass and the NN interactions are 
successful. The deviations between the exact (with charge-symmetry breaking interactions) 
and the approximate (with isospin invariant interactions) calculations are at the level of a 
fraction of a percent. The best correspondence between the results for the charge- dependent 
and charge-independent Hamiltonians is obtained for 3 H and 3 He and for all six states of 
6 Li. Moreover, we performed calculations for the T=l states in 6 Li, using the isoscalar 
potential for T=0. Analogous calculations for the T=0 states of the same nucleus were also 
performed utilizing the isoscalar potential for T=l. In all cases the absolute values of differ- 
ences between charge- dependent and charge-independent results are some three times larger 
in comparison with calculations using the correct value of the isospin for the state under 
investigation. As can be seen from the results in Table I, the biggest difference between 
energies is obtained for 6 Be, where it equals 1-2% of the "exact" value, and, to some extent, 
in 6 He, where it does not exceed 1%. 
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VI. SUMMARY AND CONCLUSIONS 



As is well-known, charge symmetry in nuclei is broken at the very least due to the 
difference between the proton and neutron masses, the charge-dependence of the strong 
interaction and the Coulomb interaction among the protons. These effects are not very 
large for light nuclei in comparison with the binding energy of the entire nucleus, but there 
are cases, in which they are significant and play a very important role, such as for nuclei 
near drip-lines, which are loosely bound. However, the description of nuclei using a fully 
charge- dependent Hamiltonian is complicated, so the construction of a Hamiltonian, as close 
as possible to charge-independent one, would be useful. 

The present article is devoted to the consideration of this problem. It is shown that 
charge-dependence of the realistic Hamiltonian for light nuclei can be taken into account 
by replacing the realistic one by an equivalent charge-independent Hamiltonian. The pa- 
rameters of such a Hamiltonian (i.e., the nucleon mass and the NN potential) are nucleus 
and state dependent. Our results can be formulated as simple recommendations for the 
construction of such an isoscalar Hamiltonian: 



1. Take the nucleon mass equal to 

1 . 

m = — (Nm n + Zrrip) . 



2. Take the realistic NN potential in states with two-nucleon isospin t = as they are 
defined from the neutron-proton data analyze. 



Take the isoscalar NN potential in states with t = 1 as a linear combination, Eq. (p2|), 
of the neutron- neutron (m t = 1), neutron-proton (m t = 0), and proton-proton (with 
Coulomb potential included) (m t = — 1) realistic potentials with coefficients c tmt , given 
in Eqs. 



The eigenfunctions of such a Hamiltonian have the isospin T of the nucleus as good a 
quantum number, but are different for the individual members of the isospin multiplet due 
to the dependence of the isoscalar Hamiltonians on Mt- 

In general, the results obtained are self-consistent and do not conflict with the simple 
requirements mentioned in the Introduction. Moreover, the isoscalar part of the Coulomb 
interaction has a nicer form than the earlier result suggested in ||16|| , Eq. (2.21), because the 



isoscalar Coulomb interaction of the nucleons in the present work equals zero in two-nucleon 
states with isospin t = 0. Consequently, the normalization of this interaction to the number 
of proton pairs occurs naturally. 

Our results can be considered to be quite good in all cases considered, because, due to 
our normalization conditions, Eq. fl33|) , we really have only two free parameters for the 
construction of the isoscalar interaction. The approximation given here is a good starting 
point for a perturbational account of nuclear charge- dependent effects not included in the 
isoscalar Hamiltonian. 
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VIII. APPENDIX A 



The orthogonality of the transformation Eq. (|14D follows directly from the observation 
that an arbitrary Jacobi tree can be written as a system of two subtrees, Fig. 3. The 
first (left) subtree with first-degree vertices rimi, r 2 m 2 , r q m q defines the Jacobi vari- 
ables £2,£3,---,6j , the second (right) subtree contains r 9+ im ?+ i, r g+2 m g+ 2, r q+t m q+t and 
£ g+ i,£ 3+ 2, ...,£ g+ t_i , correspondingly. Obviously, it is possible to introduce in the first and 
second cases the additional Jacobi variables proportional to the center-of-mass coordinates 
of the subtrees, defined as 

q ^ q+t 



i=l V ■■-H-r^-H-r- i= q+ l 



\J m \...q , = 1 ^ m q+l...q+t 



mix.;. 



The indices of these coordinates are negative, according to the rule that the arbitrary center- 
of-mass Jacobi variable situated near an intrinsic Jacobi variable £ Q is labelled as (£i 
and £oj £2 and £_i , and ). Suppose, the transformations from rimi, r 2 m 2 , v q m q to 
6-1, 6, £3, i q and from r q+1 m q+1 , r q+2 m q+2 , r q+t m q+t to £ q+1 , £ ?+2 , £ q +t-i are or- 
thogonal. Then the complete transformation is also orthogonal, because the transformation 
matrix from £_ q to £ , £1 is orthogonal: 

6 

Following this process to the top of the tree, we can arrive at orthogonal two-dimensional 
transformations. Therefore, the complete transformation (q + t) x (q + 1) is orthogonal. 




IX. APPENDIX B 

The independence of the matrix element, Eq. flTB]), (£JUM T |V? | £JIIM T ), of a, the 
intrinsic Jacobi variable number (a = 1, 2, A — 1) , can be shown as follows. By definition, 
this element equals the product of sums over spin and isospin variables and of integrals over 
intrinsic Jacobi variables 

X X J ' dl11 J dT l2--- f d VA-l%jnM T (VlV2---VA-l(TlV2---(TATl,T2...T A ) 

<T1,(T2...I7A T\,T2—TA 

x V; ljS>smM T {r}ir}2...r)A-i<ri,<72-(7ATi, t 2 ...ta) , 
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where ^sjum t (•••) is an intrinsic wave function of an atomic nucleus. This wave function 
is antisymmetric with respect to all nucleons. However, this property in Jacobi variables is 
not expandable in a simple way, because permutations of one-nucleon variables rj generate 
orthogonal transformations of Jacobi variables. We can simplify by multiplying the given 
expression by the integral 

r d Vo e* (r] ) e (770) = 1, 

where (770) is an arbitrary, normalized wave function of the center-of-mass. The new wave 
functions 

© (?7o) ^smM T {j]\f] 2 ...'r] A -\(J\Q 2 ...a A Tx, r 2 ...T A ) 

can be rewritten as usual antisymmetric functions of one-particle variables 

^ejuMt (rir 2 ...r A a 1 a 2 ...(T A Ti,T2...T A ) . 

Due to the orthogonality of the matrix a , Eq. (0), the integrals can be transformed into 
one-particle variables: 

drjo J drjt J dr) 2 ... J dr} A - X = J dr 1 J dr 2 ... J dr A . 

The operator can be written in the form 

A A A 

= 2J a a ,ia a ,j (Vi • Vj) = 2J «W«wV- + 2J a a,ia a j (Vi • V,) . 

4j'=l 4=1 i t j=i(i^j) 

The mean value of this operator is simply 

( V L> = E E J dv l J dv 2 ... J dr A §l mMT (r l r 2 ...r A o 1 a 2 ...o A T 1 ,T 2 ...T A ) 

aWI — VA T\,T2—TA 



X 



V^^jra/t (rir 2 ...rA(Ti,cr2...(JAri, r 2 ...r A ) 



.4 



4=1 ij=l(i¥=3) 

Due to the antisymmetry of the wave functions, the matrix elements (Vf ) and ((V» • V,-)) 
are the same for all possible values of z and j. Thus, 



A 



X) = (v A )j2**^i+ ((v A -i-v A )) y, 

4=1 ij=l(i^j) 

'A A A 

M> + <(Va-i • V A )> 



i=l j'=l i=l 



(Va> - <(Va-i • V A )> 



because all sums are independent of a , see Eqs. (|9|) and (|10|). Therefore, (V^ a ) is indepen- 
dent of a (a = 1, 2, A — 1). 
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X. APPENDIX C 



The relation between potentials, defined in Eq. holds also for the Reduced Hamil- 

tonians (RH) of the different two-nucleon states 

h °^t = j2c tmt h^ tmt 

mi 



and their matrix elements 



■;ut 



because the partial-wave expressions for the kinetic energy operators are the same in all 
cases. For this reason, the difference of a diagonal matrix element of the Hamiltonian H and 
an eigenvalue of the Hamiltonian H° for an arbitrary eigenstate of H° equals 



(H - H°) m ™ T = A ( A ~^ 



^ 1 

E,£',jTTtmt 



ijntmt J Qe 1 ,jntmt;E,jntmt \£>J 1 TV iV/'/') i if 'I \ 

£ ' £ ' 1 -s E/>ECtmt Q £>j7Vt (ejuT) f- 



There are nondiagonal contributions to the charge-dependent density matrix, because 
the eigenf unctions of the three RH operators h^ tmt (m t = 0, ±1) are different and are not 
equal to the eigenfunctions of the isoscalar RH /i '- 5 " 1 "*. The matrix of the isoscalar RH in this 
case is diagonal, i.e., 

mt 

The left-hand side of Eq. (CI) is the sum of traces of products of two symmetrical 
matrices, the first of these is the matrix of the RH and the second is the reduced density 
matrix. As is seen from (CI), this sum contains only two kinds of terms - one is equal to 
the product of diagonal elements of both matrices, and the other is equal to the product 
of nondiagonal elements. The last subsum is independent of the choice of the values of 
variational coefficients Ctm t - I n an y case, its contribution is negligible in comparison with 
the first subsum, which contains diagonal elements. This result is caused by the following 
reasons: The nondiagonal elements of the charge- dependent RH matrices are negligible 
in every two nucleon state, because charge-independence and charge-symmetry breaking 
forces in light atomic nuclei are significantly weaker than original charge-independent strong 
interaction between nucleons. The nondiagonal entries of charge-dependent density matrix 
are also small in absolute value, because eigenfunctions of the charge-dependent RH are 
close to eigenfunctions of charge-independent RH for the same reason as given previously. 

The most important diagonal part of the sum in Eq. (CI) is 

Ke mt {Qs,jntm t (SJUTM T ) - C tnlt Q £ , J7rt (S JUT)} . (C2) 

e,jirtmt 
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The problem is to choose the coefficients Q mt in such a way that this sum will have its minimal 
value. Let us first recall some definitions. The distributions of probabilities, given in Eqs. 
(p3|) and (|30|), differ. In the last expression they are comparable, because corresponding 
density matrices are created by the same wave-function. The relationship between these 



probabilities, following directly from the definition of the density matrix, is (see ]nj] and 
0): 

QeJ«t mt (SJUTM T ) = Q ejnt (SJUT) . (C3) 

mt 

This equality means that the quantity on the right-hand side is a marginal probability (see, 
for example, In the case t — , all entries of (C2) vanish, when coo = 1- This value 

is consistent with the definition of the isoscalar interaction, Eqs. (H), (|33"1), because, as 
mentioned earlier, in this case we have no other choice, excluding the best one. In the case 
when t = 1, the situation is not so simple, because only the sum of the three multipliers of 
hi n * mt , corresponding to the three possible values of the isospin projection m t , vanishes, due 
to relations (|33|) and (C3). In general, for any value of t the relation: 

J2{Qe,jntm t (£JUTM T ) - c tmt Q £jwt (SJUT)} = 

m t 

is satisfied. 

This relation is an identity, and, hence, is not sufficient to define the coefficients, but 
gives us more choices for determining the minimal value of the sum given in Eq. (C2). 
The best choice occurs when the partial sum of the coefficients over all possible two-nucleon 
channels is equal to zero, namely, 

{Qe^tmt {SJUTM T ) - CtmtQejirt (SJUT)} = 0. 

This equation can also be expressed in the form 

Q tmt (TM T ) = c tmt Q t (T). (C4) 
The marginal probabilities are defined by 

Y,Q^nt mt (SJUTM T ) = Q tmt (TM T ) 

and 

J^Qtmt (TM T ) = Q t (T) 

mt 

and can be extracted from the expressions for the simple multiparticle operators with known 



eigenvalues, such as the mass number A, the charge Z and the total isospin T, (see JTT 
|T3|). Therefore, 



18 



Qi CO 
Q!_! (TM T ) 
Qio (TM T ) 
Q 1+1 (TM T ) 
Using Eq. (C4), we obtain 



Qoo [1 Mt) = Qo [1 ) (A - 1) : 

3A (A - 2) + AT (T + 1) 
AA {A -I) ' 
(A - 2M T ) (A - 2M T - 2) 
AA (A - 1) ' 
A (A - 2) + AT (T + 1) - 8M| 
(A - 1) ' 
(A + 2M T ) (A + 2M T - 2) 
4A (A - 1) ' 



c tmt =Q tmt (TM T )/Q t (T). (C5) 
Applying this expression, we obtain the values of the coefficients given in Eqs. fl34p-(l37|)- 
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TABLE I. Comparison of energies, in MeV, obtained in the no-core shell-model calcu- 
lations with the full charge-dependent effective interactions, Ep_jv, and the corresponding 
isospin-invariant effective interactions, Ejso ( a s defined in the text), for A=3 and A=6 
nuclei. The effective interactions were derived from the CD-Bonn potential ||J]. Harmonic- 
oscillator energies of HQ =19 MeV and Ml =17.2 MeV were employed for the A=3 and A=6 
systems, respectively. For A=3 a complete 8hfl model-space was utilized, while for A=6 a 
complete 6HQ model-space was used. 



A Z(J n T) 


Ep_Ar 


E/so 


[Ep-N — Eiso] / Ep_N 




-8.441 


-8.421 


-2.44 x 10~ 3 


3 He(i + |) 


-7.668 


-7.643 


-3.27 x 10~ 3 


6 He(0 + l) 


-25.665 


-25.830 


+6.42 x 10~ 3 


6 He(2+l) 
6 Li(l+0) 


-22.944 


-23.085 


+6.13 x 10~ 3 


-28.257 


-28.180 


-2.73 x 10~ 3 


6 Li(3+0) 


-25.605 


-25.525 


-3.12 x 10~ 3 


6 Li(0+l) 


-24.815 


-24.666 


-6.02 x 10~ 3 


6 Li(2+0) 
6 Li(2+l) 
6 Li(l+0) 


-23.146 


-23.071 


-3.25 x 10~ 3 


-21.934 


-21.924 


-4.47 x 10~ 4 


-20.433 


-20.362 


-3.46 x 10~ 3 


6 Be(0 + l) 
6 Be(2+l) 


-22.829 


-22.545 


-1.25 x 10~ 2 


-20.227 


-19.851 


-1.86 x ID" 2 
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A. Figure captions 

FIG. 1. An example of a Jacobi tree. 

FIG. 2. The Jacobi tree for a nucleus composed of Z protons and N neutrons. 
FIG. 3. An example of a Jacobi tree with subtrees. 



22 




FIG. 1 . The Jacobi tree sample. 

Isoscalar Hamiltonians for light atomic nuclei 
G.P.Kamuntavicius, B.R.Barrett, P.Navratil, 
G.Sapragonaite, and R.K.Kalinauskas 




FIG. 2. The Jacobi tree for nucleus composed of Z protons 
and N neutrons. 
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FIG. 3. The subtrees of Jacobi tree. 
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